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Using exact diagonalizations, Green's function Monte Carlo simulations and high-order perturba- 
tion theory, we study the low-energy properties of the two-dimensional spin-1/2 compass model on 
the square lattice defined by the Hamiltonian H = — ( J^^ir^cr^+e^ -I- Jzcr^ir^+e^ ). When J^c Jz, 
we show that, on clusters of dimension L x L, the low-energy spectrum consists of 2^ states which 
collapse onto each other exponentially fast with L, a conclusion that remains true arbitrarily close 
to Jx = Jz- At that point, we show that an even larger number of states collapse exponentially 
fast with L onto the ground state, and we present numerical evidence that this number is precisely 
2x2^. We also extend the symmetry analysis of the model to arbitrary spins and show that the 
two-fold degeneracy of all eigenstates remains true for arbitrary half-integer spins but does not apply 
to integer spins, in which cases eigenstates are generically non degenerate, a result confirmed by 
exact diagonalizations in the spin-1 case. Implications for Mott insulators and Josephson junction 
arrays are briefiy discussed. 

PACS numbers: 71.27.+a,75.30.Ds,03.67.Lx 



I. INTRODUCTION 

Building on the deep understanding of the Heisenberg 
and other models of magnetism, it is a very common prac- 
tice to describe discrete degrees of freedom as pseudo- 
spins, with the hope to gain insight from the form of 
the resulting magnetic model. A well-known example of 
considerable current interest shows up in the context of 
Mott insulators with orbital degeneracy: In an octahe- 
dral environment, the degeneracy of the d electrons is 
only partially lifted, and the remaining orbital degree of 
freedom is often described as a spin-1/2 or 1 for eg and 
t2g electrons, respectively However, the situation 

is in general less simple that one might hope. Indeed, as 
already emphasized by Kugel and Khomskii, the symme- 
try of pseudo-spin Hamiltonians is in general much lower 
than SU(2), and there are cases where the properties of 
the resulting model are poorly understood. This is in 
particular the case of models where the anisotropy of the 
coupling in spin space is related to the orientation of the 
bond in real space. The simplest version of such a model 
on the square lattice is defined by the Hamiltonian 



r r+ej; 



(1) 



where r" are the x and z components of a pseudo-spin 
operator. By analogy with the dipolar coupling between 
compass needles, this model has been called the compass 
model by Kugel and Khomskii 0. Realistic models of 
orbital degeneracy are usually more complicated in sev- 
eral respects. In particular, the spins and pseudo-spins 
are usually coupled. Pure orbital models can be of direct 
relevance though if the spins order ferromagnetically, as 
recently argued by Mostovoy and Khomskii in the con- 
text of NaNi02 [3 . The precise symmetry is also usually 
more complicated than the simple case of this Hamilto- 



nian, but we will nevertheless concentrate on that model, 
considering it as a minimal model rather than a realistic 
one. 

Interestingly, such models have appeared in other con- 
texts as well. First, even in Mott insulators without or- 
bital degeneracy, extra degrees of freedom can appear 
if the system is frustrated and, for instance, consists of 
spin-1/2 coupled triangles, like in the trimerized kagome 
lattice, in which case the chirality that keeps track of 
the extra degeneracy of each triangle plays a role similar 
to that of orbitals In a magnetic field, this model 
has been predicted to exhibit a magnetization plateau at 
1/3 1^, and the low-energy properties inside the plateau 
can be described by a kind of compass model . 

More recently, the model of Eq. with anisotropic 
couplings along the x and z directions has been proposed 
by Dougot and collaborators in the context of Josephson 
junction arrays 0. 

Despite its deceptive simplicity, the model of Eq. 
is a formidable challenge, in many respect comparable 
to very frustrated magnets. To see this, let us follow 
Ref. Q and consider the classical version of the model, 
in which spins are considered as classical vectors. In that 
shown by Nussinov et ai, the ground state is 
highly degenerate, as in very frustrated magnets. First 
of all, all ferromagnetic states are degenerate, regardless 
of the relative orientation of the spins with respect to the 
lattice, as can be easily checked from Eq. . In addition, 
from any ferromagnetic state, one can construct other 
states by flipping all spins of a z column with respect to 
an X mirror, or equivalently by flipping all spins of an x 
line with respect to a z mirror. Since all these operations 
can be performed simultaneously and in any order, they 
generate a discrete degeneracy of order 2^. 

The effects of thermal fluctuations on the classical 
model have been convincingly identified by analytical and 
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numerical approaches. Nussinov et al. have shown that 
an order by disorder mechanism is expected to hft the ro- 
tational degeneracy and to select states in which the spins 
point along the x or z axis, leading to a nematic ground 
state since lines or columns of spins are still free to flip. 
Using extensive Monte Carlo simulations, Mishra et al. 
have shown that the two possible orientations along x or 
z lead to an effective Ising order parameter, and that the 
model undergoes a finite temperature phase transition of 
the Ising type . 

On the other hand, the understanding of the spin-1/2 
version of the model is still preliminary. Most of the 
results have been obtained by Dougot et al. in their 
analysis of a generalized version of the model of Eq. ^ 
defined by the Hamiltonian 

i? = - E ( J. <t;+,^ + J, ) , (2) 

r 

in which the couplings along the x and z directions 
can take different values. Using elegant symmetry argu- 
ments, DouQot et al. have shown that all eigenstates must 
be at least two-fold degenerate. They have also shown 
that in the strongly asymmetric case {Jx/ Jz ^ 1 or 
Jz/Jx ^ 1), the 2^ states that evolve adiabatically from 
the 2^ ground states of the decoupled Ising chain case 
{Jx — or Jz — 0) should collapse onto each other. This 
perturbative argument does not apply close to Jx = Jz 
though, and whether this remains true in the isotropic 
case could not be decided. Note that this is an impor- 
tant issue in the context of quantum bits in which they 
came accross this model since the presence of a gap would 
help to protect the q-bits. More recently, Nussinov and 
Fradkin ^3 have shown that these models are dual to 
models of p + ip superconducting arrays, and have dis- 
cussed general properties of order parameters and phase 
transitions. 

In this paper, we concentrate on the zero-temperature 
properties of the quantum version of the model. We 
start by a semi-classical analysis of the model of Eq. 
and show that, as in many frustrated magnets, quantum 
fluctuations essentially have the same effect as thermal 
fluctuations regarding the lifting of classical degeneracy 
(section!^. We then turn to an extensive analysis of the 
model of Eq. ^ in the spin-1/2 case. As we shall see, 
there is no minus sign problem in the Green's function im- 
plementation of quantum Monte Carlo, which allows one 
to study the ground state properties of very large clus- 
ters, i.e., up to 17 X 17. Combined with the large number 
of symmetries, hence of different symmetry sectors that 
can be studied independently, this allowed us to reach 
definite conclusions regarding the low energy spectrum, 
conclusions that agree with high order perturbation the- 
ory. These results are presented in section Hill Finally, 
the symmetry analysis of Dougot et al. is extended to 
larger spins in section llVl with the conclusion that integer 
and half-integer spins behave once more quite differently. 
Some implications of the present results are discussed in 
the last section of the paper. 



II. SEMI-CLASSICAL COMPASS MODEL 

In order to have a first insight into the properties of 
the quantum version of the compass model, we have per- 
formed a spin-wave analysis in the symmetric case de- 
fined by Eq. In that respect, it is useful to emphasize 
that, as noticed in Refs. [1,0, the degeneracy is partly 
accidental and partly due to symmetry. Indeed, in addi- 
tion to the lattice translational symmetries, this model 
has two types of discrete symmetries: (i) The Qi trans- 
formation which flips the z component of all the spins 
of the column — i, and the Pj transformations which 
flip the X component of all spins of the line = j. Note 
that the transformations Qi (resp. Pj) could be seen as 
rotations of all the spins of the column (resp. line) about 
the Bx axis (resp. e^) by an angle tt. (ii) The simultane- 
ous rotation i?y(f ) of all spins and of the lattice about 
an Ey axis by an angle ^. Starting from any state, these 
symmetries generate new states with exactly the same 
classical energy and the degeneracy associated to these 
symmetries cannot be lifted by thermal fluctuations 
However, the ground state also has an accidental degen- 
eracy: All ferromagnetic states are degenerate regardless 
of the angle between the spins and the lattice. This de- 
generacy is not related to a symmetry since the model is 
not rotationally invariant. Accordingly, it has been found 
that the thermal fluctuations partially lift this degener- 
acy via an order-by-disorder mechanism, favouring the 
2x2^ ground states with all spins parallel to Bx or 
and implying a directional ordering of the spins. These 
favoured states are the uniform state with all spins in 
the direction Bx and all the states obtained by applying 
the symmetries Ry(^), Pj and Qi. So, as anticipated, 
only the accidental degeneracy is lifted by thermal fluc- 
tuations. 

The same idea applies to quantum fluctuations. Start- 
ing from an arbitrary ground state, one can bring it back 
into a ferromagnetic ground state applying only sym- 
metry operations. But applying symmetry operations 
does not change the form of the Hamiltonian. The fluc- 
tuations around both states will thus have exactly the 
same form. Then it is sufficient to do the spin-wave 
expansion around the uniform classical ground states 
Sr = ScosOqBx + SsuiOqBz. To linear order in 1/S, 
the energy can be brought into the form 

Q 

with LUq{6Q) — 4 JS'V 1 — cos^ 6*0 cos kx — sin^ 0o cos k^, 
and the resulting ground-state energy is plotted in Fig. Q 
As in the classical case with thermal fluctuations, the fig- 
ure clearly shows that the angles = 0, §, tt, ^ are se- 
lected by quantum fluctuations since they minimize the 
energy. Applying the symmetries Qi, Pj and -Ry(-|) to 
these states gives 2x2^ equivalent favoured states, corre- 
sponding to 2^ states parallel to Bx and 2^ states parallel 
to Bz. So, as in the classical case, there is a directional 
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FIG. 1: Ground-state energy of the ferromagnetic states de- 
fined by Sr = S{9o) = ScosOoBx + SsinOoBz including zero- 
point energy as a function of the angle ^o- Ed = ~JNS^ is 
the classical ground-state energy. 

ordering of the ground state. 

As is many frustrated magnets like the J1 — J2 model on 
the square lattice for J2/J1 = 1/2, this calculation is not 
fully consistent since the correction to the magnetization 
diverges the present case, the divergence comes 

from a line of zero energy along the = direction when 
00 = 0, TT and along the fcz = direction when = t"- 
We have pushed the expansion to next order in 1/S and 
checked that a self-consistent mean-field decoupling of 
the 4-boson terms suppresses the divergence. At this ap- 
proximation, the spectrum becomes gapped, which does 
not violate any general theorem since the model does 
not possess rotational symmetry, and long-range order is 
preserved even for spin-1/2. 

In the anisotropic case {Jx ^ Jz), discrete rotational 
symmetry of the Hamiltonian as well as continuous ro- 
tational symmetry of the ground state are lost, leading 
to only the 2^ ground states with rows of parallel spins 
along the x axis for > Jz and columns of parallel spins 
along the z axis for < Jz- Since these 2^ classical 
ground states are related by symmetries of the Hamilto- 
nian, quantum fluctuations cannot lift this degeneracy. 



III. QUANTUM COMPASS MODEL: SPIN-1/2 

In this section, we turn to the spin-1/2 case for which 
we write the Hamiltonian 

H = -Y, {Jx«+.^ + J^KK+a^) : (4) 

r 

where and ct^ are Pauli matrices acting on the spin at 
site r. In the following, we will use the parametrization 
of the exchange integrals Jx = J cos 9 and Jz = J sin 9 
with 9 G [0, and we will study the model on A^-site 
square clusters of dimension L x L. 



Taking all symmetries into account, exact diagonaliza- 
tions could be performed up to L = 5. As we shall see, 
this is not sufficient to draw conclusions regarding the de- 
generacy of the ground state in the thermodynamic limit 
close to Jx = Jz- However, this model has the very inter- 
esting property that all non-diagonal matrix elements are 
negative. This has allowed us to implement the Green's 
function Monte Carlo algorithm 12], which gives access 
to the ground-state energy in a given symmetry sector, 
and to reach clusters up to L = 17. Besides, one can 
choose quantum numbers so that all relevant low-energy 
states are ground states of a given symmetry sector. To 
see how this works, let us look more closely at the sym- 
metries of the model. 

In addition to the lattice translation symmetries, the 
Hamiltonian |0J has another type of discrete symmetries. 
The first one corresponds to the operators Qi — o'fj i 
which are the products of the tr^ on one column (tx — i)- 
These operations correspond to a rotation by an angle 
TT about the e^, axis of all the spins of a given column: 
Q'^afjQi = -erf/ and Q^^a-f jQ^ = af^ The second 
one corresponds to the operators Pj = Yii '^i which are 
the products of the on one line (r^ — j), and which 
correspond to a rotation by an angle tt about the axis 
of all the spins of a given line: P^^a^'j'Pj — —cr^f and 
Pj'^cr^jPj = crf^ . In the isotropic limit Jx = Jz, there is 
one more discrete symmetry: the global rotation Ry{^) 
of all spins and lattice about the Sy axis by an angle ^. 

As emphasized by DouQot et al, the Pj's commute 
with each other, as well as the Q^'s, but [Qi,Pj] ^ 
y This has two remarkable consequences: First of all, 
all eigenstates must be two-fold degenerate. Besides, we 
can choose either the Pj 's or the Qi's to define symmetry 
sectors in which the Hamiltonian can be independently 
diagonalized. Since Pj — 1 and = Pj, the eigenvalues 
of Pj are pj = ±1, the same being true for the Qi's. Thus 
the Hamiltonian can be diagonalized in the symmetry 
sectors characterized by the set (pi, • • • Pi = ±1, or 

alternatively in the sectors defined by the eigenvalues of 
Qi and characterized by the set (gi, • • • , g^), gi = ±1. 

To see how this works, let us start from the trivial case 
Jx = 0. The model then consists of a set of decoupled 
Ising columns with eigenstates jm) = |mi^i) ® lmi.2) (8) 
•■•(g) Im^.L), where mi,j — ±1 is the eigenvalue of af j. 
The ground state manifold contains 2^ states defined by 
rui^i — -- - = TTii^L = il, i = 1,---,L. Clearly, all |m) 
states are eigenstates of the Pj's- Now, in any of the 
ground state, all pj's are equal since all lines are identical, 
and the ground state manifold consists of 2^~^ states in 
the sector pi ^ ■ ■ ■ = pl = +1 and 2^~^ states in the 
sector pi = •■• = PL = — 1^ 

The classification according to the Qi's is quite differ- 
ent. First, note that the states {|»ti)} are not eigenstates 
of Qi. However, if we denote by | t)i (resp. | l)i) the 
ground state of column i with all the spins up (resp. 
down), then one can define two new ground states by 
= 45(IT), + U)0 and = 1 (It)i-U),)^ 
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These new states are eignestates of Qi with eigenvalues 
±1 respectively. So qi can be either —1 or +1 for each 
column independently, and the ground state manifold 
has one member in each of the 2^ sectors (gi, • • • ,(7l), 
Qi = ±1. 

The case Jz = is connected to the case Jj. = by 
the rotation symmetry Since Ry{^)^^ PjRy{^) — 

{-1)^Ql-3 and Ry{%)-^Q^Ry{^) = the symme- 
try sectors are interchanged. So the ground state is 
made up of 2^ states in each 2^ sectors (pi, • • ■ = 
(±1,... ,±1). 

So, as announced earlier, to determine the structure of 
the low-energy spectrum when going away from — Q 
or = 0, it is always possible to choose the quantum 
numbers so that each state is the ground state of a given 
symmetry sector. 



A. Exact diagonalization 

The spectrum of the Hamiltonian versus the asym- 
metry parameter 9 has been determined for L = 2,3,4 
and 5 using exact diagonalization. Fig. |21 presents the 
low energy levels for L = 4 and L = 5. As expected, 
when Jx,Jz 7^ 0, the ground state is two- fold degenerate 
with one state in the sector Pi = -- -= pl = +1 and the 
other in the sector pi = ■ ■ ■ = pj^ = —1. For < < J^, 
the 2^-fold degeneracy of the ground state at 6* = (i.e., 
Jz ~ 0) is lifted by the Jz term of the Hamiltonian. The 
lowest of these states is in the sector pi = ■ ■ ■ = pi^ = ±1 
and the highest in the sector pi = ±( — 1)*. The gap be- 
tween these two states is denoted by A(0). As we will 
see in the following, A(6') goes to zero for all in the 
thermodynamic limit. 

In the symmetric case {9 = i.e., Jx = Jz), the 
adiabatic continuation of the degenerate ground states 
of the Jx = and = cases generates only 2 x 2^ — 2 
states, and not 2 x 2^ as one might naively expect from 
the semi-classical case. The reason is that the lowest pair 
of state is common to the two families of states coming 
from Jx = and Jz — 0, while all other states cross at 
Jx ~ Jz (see Fig. 12). This does not mean however that 
the low energy sector has only 2x2-^ — 2 states for very 
large systems. In fact, another pair of states is decreasing 
quite fast toward the ground state as a function of the 
size. Since these states correspond to the first two-fold 
degenerate excited state of the sectors pi = • • • = _pl = 1 
and pi — ■ ■ ■ ^ PL ~ —1, one can keep track of their 
energy for all sizes. The gap between the ground state 
and these states, denoted by A2, is plotted in Fig. O as a 
function of 1/N. These results are indeed consistent with 
a vanishing of this gap in the thermodynamic limit [l3j | . 



B. Perturbation theory 

Before discussing the results obtained for large clusters 
with Green's function quantum Monte Carlo, let us see 
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FIG. 2: Energy of the low- lying states versus anisotropy pa- 
rameter 6 for 4x4 and 5x5 lattices obtained by exact diag- 
onalization {Jx = J COS 61, Jz = J sin 6'). 



what perturbation theory predicts regarding the scaling 
of the gap A in the limit Jz Jx- To this purpose, the 
Hamiltonian is written as 



H = Ho + V, 



(5) 



where Hq = - Jx J2r ^r'^r+e^ ^ud V = - Jz J2r '^r^r+e, ■ 

Since each term of V flips a pair of spins on neighboring 
chains, it is necessary to apply the perturbation at least 
L times on an L x L lattice to flip all spins of a pair of 
neighboring chains and reach another ground state. So 
the gap is expected to scale as A/Jx = a{bJz/Jx)^ ■ One 
can be more precise though and determine the constants 
a and b from high order perturbation theory. To first 
significant order in Jz/Jx, the gap A is given by the 
gap between lower and higher eigenvalues of the effective 
Hamiltonian 



P^VP^.V ■■■VP,,_,VP, 



Ui>0 VL-1>0 



(6) 

where P^ is the projector on the v-th eigenspace of Hq. 
This eigenspace has an energy E^, = —JxN + AvJx, where 
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FIG. 3: Gap A2 (see text) as a function of the number of sites 
iV. 



nothing but the Hamiltonian of the one-dimensional Ising 
model. So, the gap between the lowest and highest 
eigenvalues is Amax — ^min — 2L when L is even and 
Amax — Amin — 2(i — 1) wheu L is odd. Finally, if we 
define P{L) by 



1 



(7) 



the gap A becomes 



%LP{L)(J^/AJ^)^ if L is even 

A/ J, = { (8) 
8(L - l)P{L){J^/A.J,^)^ if L is odd. 

The dominant behaviour of P{L) has been determined 
numerically (see Fig. 0)). It turns out that LP{L) ~ 
exp(0.754i - 0.694), which leads to: 



1/ can take the values v — 0,1, ■ ■ ■ , L/2 when L is even 
and 1/ ~ 0,1, ■ ■ ■ , {L — l)/2 when L is odd. Let {|;/, k) : 
fc = 0, 1, • • • } be a basis of the i^-th eigenspace of i?o- We 
have to evaluate 



(0,l\VP,,V---VP,^_,V\0,k) 



■p. 



(0,- 



z z r> 



For L > 2, the only contribution to A arise when the 
a product flips two neighboring lines with = j and 
= j + 1. It corresponds to ji = • • • = Jl = j and 
ii — p{l), ■ ■ ■ ,iL = p{L) with p ^ Sl and Sl the set of 
L\ permutations. The ground state |0, k) contains only 
ferromagnetic lines in the x or —x direction. So, for m < 



L, the state ct^, , , , ■ , , ■ 



•'^p(L)jVi)j + 



ilo,fc) 

cannot be in the ground state manifold of H^, but it must 
be in one of the excited eigenspaces, say the z^m(?')-th. 
Then we can write: 

{0,l\VP,,V ■ ■ ■VP,^_^V\0,k) 



L-l 



•^p(L)j^p(L)j+iio>fc) n 



Using ~ Eq + AvJx and rearranging the cr^'s, we get: 



{0,l\H[^\0,k) 



W'-\^^^M---^L-i{p) 
(0, ;| E • • • j^i j+i |o> k). 
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FIG. 4: log-linear plot of LPiV) versus the linear lattice size 
L. The squares are the results obtained numerically, and the 
dotted line is an exponential fit. 



A/J, = 



3.997 (0.531JjJ:r) 



if L even 



3.997(1 - IjL) {OmiJ,JJ^Y if L odd. 



(9) 

It has the correct [Jz/Jx) behavior which corresponds 
to the exact result for the two- line system 0. So, when 
{Jz/Jx) 1, this approximation predicts A — > in the 
thermodynamic limit. Moreover, since Q.531{Jz/ Jx) < 1 
even for Jz/Jx = 1, it seems likely that this scaling will 
remain true up to the symmetric limit J^ — J^. As we 
shall see, this is confirmed by the Green's function Monte 
Carlo results. 



Green's function Monte Carlo 



Next, we note that, in terms of the pseudo-spin rj = 
<i • • • the operator Y.^ ' ' ' o-LjO-Lj+i is 



If the Hamiltonian of a model has only non-povitive 
off-diagonal matrix elements, which is the case here, the 
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Green's function Monte Carlo method [IJ allows one to 
calculate the ground-state energy of a given symmetry 
sector by using a stochastic approach. The algorithm we 
have used is the implementation with a fixed number of 
walkers described in detail by Calandra and Sorella [l^ . 
and the guiding function is given by: 



exp 




(10) 



where \Fx) is the ferromagnetic state in x direction such 
that <Jr\Fx) — \Fx) Vr. The parameters Vr,r' were de- 
termined in order to minimize the energy of the guiding 
function. The ground state energy has been determined 
in each sector (pi, • • • ,Pl) separately, which corresponds 
to the energy of the first 2^ states of the full Hilbert 
space (the states coming from the limit Jz — ground 
state in Fig. I^J, and gives access to the gap A(0). 

Let us first discuss the scaling of the energy per site as 
a function of the system size. As can be seen in Fig. El the 
energy per site is strongly size dependent up to a certain 
size (8x8 for = Jz), and is very little size dependent 
for larger clusters. This indicates that strong finite-size 
effects are to be expected, especially close to the sym- 
metric point, justifying the use of quantum Monte Carlo 
to get information on large clusters. 
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FIG. 6: Gap A versus linear lattice size L for various values 
of the asymmetry parameter 9. For L = 2, 3, 4, 5 the results 
were obtained from exact diagonalization. 
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FIG. 7; Q versus Jz/Jx, where a is the gap scaling constant 
A ~ Q^. The symbols are the results obtained with Green's 
function Monte Carlo and the dotted line is fitted to these 
values. 



FIG. 5: Ground-state energy Eo versus the number of sites 
A*' for various values of the asymmetry parameter 6. 

Fig. shows a log-linear plot of the gap A versus the 
linear lattice size L for various values of the asymme- 
try parameter 9. Perturbation theory predicts that the 
scaling of A is given by a power law A/J^ oc a^, with 
a ~ 0.531Jz/ Jx- To check this prediction, we have fitted 
the results of Fig. with a straight line for each value 
of 9, keeping only sizes beyond which the scaling is ap- 
proximately linear. The values of a deduced from this fit 
are plotted in Fig.[7|as a function of Jz/Jx- Remarkably, 
the relation between a and Jz/Jx is quite linear up to 
Jx = Jz, which seems to justify the perturbation the- 
ory in this limit. Moreover, a fit gives a ~ 0.55Jz/ Jx, 



in good agreement with the perturbation theory predic- 
tion a ~ 0.531 Jz/ Jx- In fact, for small Jz/Jx, even the 
prefactor and the even-odd effect predicted by the per- 
turbation theory agree with the results of Fig. All 
these results lead us to the conclusion that the gap A 
indeed follows a power law A/J ^ a^, with a < 1 when 
Jz < Jx, implying that A — > in the thermodynamic 
limit. Thus the 2^ states coming from the limit 9 = 
{Jz = 0) ground state collapse in the thermodynamic 
limit, as long as Jz < Jx- Using the rotation symme- 
try Ry{^), the same result holds true for the 2^ states 
coming from the limit 9 ~ ^ (he., Jx = 0) ground state, 
as long as Jz > Jx- The conclusion is that there are 2^ 
states collapsing exponentially fast onto each other when 
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Jx ^ Jz: and at least a 2 x 2^ — 2 when (i.e., Jx = Jz)- As 
we argued above, the actual number is very probably ac- 
tually equal to 2 x 2^, in agreement with the semiclassical 
analysis. 

IV. GENERAL SPIN 

The symmetry arguments used for spin-1/2 can be eas- 
ily extended to arbitrary spins. Let us consider the sys- 
tem with N spins S* on a L x L lattice described by the 
compass model Hamiltonian 

H — ~ (^JxS^S^^^^ + JzS^S^^^^) . (11) 

r 

It is straightforward to check that the generalizations 
of the Pj and Qi defined by 

Qi = n*e— (13) 

for integer spins and by 

Pj = ]Je-'"^^^- (14) 

3 

for half-integer spins commute with the Hamiltonian. 
However, it is easy to check that the commutator [Qi, Pj] 
vanishes for integer spins, whereas it does not van- 
ish for half-integer spins {{Qi,Pj} = 0). Indeed, the 
only terms in the product which do not trivialy com- 
mute are those on the site So we just have to 
show that [e~^'^^^-i ,e~™^'-^] = for integer spins and 
|g-i7r5i J g-47rSj^ I _ Q ];ialf-integer spins. The opera- 
tor e~"'^*'j = Rz{tt) corresponds to the spin rotation by 
an angle tt about the z axis acting at site Applying 
this rotation to e~'^'^''J gives 

i?;i(7r)e— ^?..i?,(7r) = ^ (^i?;i(^)(5-^.)"i?.('r) 



But for integer spins, e^^^^-i = e"*'^'^'-^ , while for half- 
integer spins, e^'^'^'-J = — e~"'^''J , which terminates the 
proof. This has two important consequences. First of all, 
in DouQot et al.'s argument, the fact that [Qi, Pj] ^ was 
crucial to show that each state was two-fold degenerate. 
We thus expect that for all integer spins, the eigenstates 
are generically non-degenerate. Besides, for integer spins, 
since all Qi's commute with all Pj 's, one can use all these 



symmetries simultaneously, leading to 2^^ different sym- 
metry sectors. 

We have checked these predictions for 5 = 1 with ex- 
act diagonalizations of clusters of size L — 2 and L — 3 
(see Fig. (SI . Indeed, the ground state is non-degenerate 




FIG. 8: Spectrum versus asymmetry parameter 9 for linear 
lattice sizes L = 2 and L — 3 obtained by exact diagonaliza- 
tion. Jx — J cos 9 and — J sin 9. 

for all Jx, Jz =/= 0, and it is in the only symmetry sector 
Pi ^ ■ ■ ■ ^ PL = qi = ■ ■ ■ = Ql = (-1)^ connecting 
the ground states in the limits Jx = and Jz = 0. The 
Jz = (resp. Jx = 0) ground state 2^-fold degeneracy is 
partially lifted by the Jz (resp. Jx) term of the Hamil- 
tonian, creating a gap A. As in the spin-1/2 model, the 
gap between the ground state and the first excited state 
which has the same degeneracy and is in the same sym- 
metry sector as the ground state is denoted by A2. The 
gaps A and A2 are smaller for L = 3 than for L — 2, and 
one can conjecture that these gaps go to zero in the ther- 
modynamic limit. In this case, the ground state would 
have the same degeneracy as for spin-1/2. A definite 
conclusion would clearly require to study larger clusters 
though. 
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V. CONCLUSION 

Using a variety of approaches, we have obtained a co- 
herent picture of the zero-temperature properties of the 
quantum compass model. On a finite cluster, we have 
confirmed that all eigenstates of the spin-1/2 model are 
at least two-fold degenerate, a result that we have ex- 
tended to arbitrary half-integer spins, while they are not 
necessarily degenerate for integer spins. However, the 
degeneracy that remains when thermal or semi-classical 
quantum fluctuations are introduced, namely the pos- 
sibility to flip the spins along lines or columns, is still 
present as a manifold of states which collapse exponen- 
tially fast onto the ground state upon increasing the size 
of the lattice. This was already known to be the case for 
the asymmetric case not too close to = Jz- Thanks 
to extensive quantum Monte Carlo simulations, we have 
shown that this remains true up to the symmetric case 
Jx = Jz, and that the the number of these states (2^ 
when Ja; 7^ Jz, 2 X 2-^ when = Jz) agrees with the 
degeneracy predicted by the semi-classical analysis. 

Physically, this has two consequences. Regarding or- 
bital fluctuations in Mott insulators, in which case the 
symmetric version of the model seems more appropri- 
ate, the present results confirm the absence of true or- 
bital long-range order for quantum spins in the thermo- 
dynamic limit even at zero temperature. Nevertheless, 



as in the case of thermal fluctuations for classical spins, 
the possibility to choose between the x and z directions 
should still lead to a finite-temperature Ising transition. 

Regarding Josephson junction arrays, one of the im- 
portant issues is to ensure that the two-fold degenerate 
ground state is well protected by a gap to all excited 
states. As noticed by Dougot et ai, this requires to 
work with not too large systems if a family of states 
collapse onto the ground state in the thermodynamic 
limit, as they already showed to be the case for J^ Jz- 
Their results suggested however that there might 
be a quantum phase transition to a gapped phase 
around Jx = Jz in which the gap to all excited states 
would remain finite even in the thermodynamic limit. 
This possibility is clearly ruled out by the present results. 
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